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We present a combined experimental and computational methodology for the discovery of new 
materials. Density functional theory (DFT) formation energy calculations allow us to predict the 
stability of various hypothetical structures. We demonstrate this approach by computationally 
predicting the Ce-Ir-In ternary phase diagram. We predict previously-unknown compounds Celr4ln 
and Ce2lr2ln to be stable. Subsequently, we successfully synthesize Celr4ln and characterize it by 
X-ray diffraction. Magnetization and heat capacity measurements of Celr4ln are reported. The 
correct prediction and discovery of Celr4ln validates this approach for discovering new materials. 

PACS numbers: 61.66.Fn, 71.15.Mb, 71.15.Nc, 61.05.C- 



I. INTRODUCTION 

The holy grail of computational materials design is the 
ab initio discovery of new materials which can be syn- 
thesized to have specific desired properties. Our increas- 
ing understanding of the structure/properties relation- 
ship in many materials and our increasing computational 
resources make this an exciting possibility.'^ However, in 
many cases, the precise relationship between structure 
and properties is not yet well understood. This is the 
case for strongly correlated electron materials which ex- 
hibit interesting phenomena such as superconductivity, 
colossal magnetoresistance, and multiferroic behavior.lSl 

Despite lacking a theory to immediately predict the 
properties of strongly correlated materials from struc- 
ture alone, phenomenological observations often suggest 
families of composition and types of structures which are 
likely to exhibit specific functionalities. Having identified 
a family of such structures, the question becomes whether 
one can synthesize new materials in these families having 
enhanced properties. Unfortunately, the number of can- 
didate materials which one might identify is often very 
large, and the experimental costs in searching for com- 
pounds can thus be heavy. Therefore, a computational 
method for predicting the stability of hypothetical com- 
pounds is highly desirable. 

A variety of computational structure search methods 
have been presented including random search methods, 
genetic algorithms and other heuristic methods, and so- 
phisticated data-mining approaches.l^H^l However, such 
methods are non-ideal for our application; because we 
are interested in searching for structures at many differ- 
ent compositions, random search and heuristic methods 
often require a prohibitively large amount of computa- 
tion time. Formal data-mining approaches often fail in 
ternary and higher-order systems as there is an insuffi- 
cient amount of data available. In addition, all of these 



methods require a significant amount of resources de- 
voted to software development. In our case, the parallel 
progress of experiment and theory led to fast results. Nei- 
ther theory nor experiment alone could have progressed 
to explore this material system with modest devoted re- 
sources and on the same timescale of several months. 

In this paper we report on our exploration of the Ce- 
Ir-In ternary phase diagram. Cerium-based transition 
metal indide compounds adopt a large variety of struc- 
ture and bonding arrangements which produce a large 
variety of interesting physical phenomena including su- 
perconductivity.'' Two types of structures which we are 
interested in exploring include systems with Ce sublat- 
tices which exist in either square planar or geometri- 
cally frustrated arrangements. Experience with cuprates, 
iron-pnictide, organic, and heavy fermion superconduc- 
tors suggests that the square-planar compounds are likely 
to exhibit superconductivity.^ In geometrically frustrated 
systems, large degeneracies can lead to novel physical 
phenomena.'Sl 

Using this method, we examine the structural stability 
of selected compounds in the Ce-Ir-In ternary phase di- 
agram and predict two previously unknown compounds, 
Celr4ln and Ce2lr2ln, to be stable. Subsequently, we 
were able to synthesize Celr4ln, which possesses a geo- 
metrically frustrated Ce sublattice. 

II. METHODOLOGY 

A. Thermodynamics 

At thermal equilibrium, a compound is stable if it has a 
negative formation energy with respect to its constituent 
elemental and any other competing phases. The forma- 
tion energy per atom for a ternary compound Ce„lrmln; 
with respect to its constituent elements may be written 
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FIG. 1: (color online) Predicted Ce-Ir-In phase diagram. The 
top-down view of the convex hull of the energy of the various 
Ce-Ir-In phases illustrates which phases are stable and un- 
stable at zero temperature. The ternary phases at corners of 
triangle represent stable ternary compounds, i.e. compounds 
such as Ce5lr2ln4 (1), Celrlns (3), and Celr4ln (5) are pre- 
dicted to be stable. Compounds such as Celrln2 (7), Celr2ln7 
(12), and Celr5ln2 (13) fall within the hull, rather than on its 
surface, and are predicted to be unstable. The binary and 
ternary compounds are numbered following their order in Ta- 
ble m 
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where Eqc, Ei^, and Ein are the Gibbs free energies per 
atom of the elemental phases. 

Using these formation energies, stable and competing 
phases can easily be illustrated. As shown in Fig.[l] each 
compound is plotted according to its composition and 
formation energy. Then, the smallest surface encompass- 
ing all points is known as the convex hull, and points 
which lie on the convex hull itself are thermodynamically 
stable. 

Both our computational and experimental techniques 
have limitations in finding these low-energy structures. 
Identifying candidate structures is an open challenge to 
the computational methods which also have known sys- 
tematic errors. The annealing process in experimental 
synthesis techniques may fail to find the lowest energy 
structure due to kinetic barriers and insufficient driving 
forces. Thus, synthesis routes may exist which allow us 
to produce metastable compounds and at the same time 
it may be difficult to realize the true ground state. 



B. Details of the Calculations 

Calculation of the entire Gibbs free energy is computa- 
tionally expensive and unnecessary. The internal energy 
contribution to the free energy dominates the enthalpic 
and entropic terms, particularly at low temperature and 
pressure. In this work, energies are calculated and evalu- 
ated in the zero temperature and pressure limit, neglect- 
ing zero point energies which are small in these heavy 
elements. 

Our DFT calculations are performed using 
with the projector-augmented wave PAW method in the 
generalized gradient approximation PBE.^'^ Energy dif- 
ferences between compounds were converged to within 
0.3 meV/atom using a plane wave cutoff energy of 400 eV 
and a k-point mesh of density 40 A~^. During the struc- 
tural optimization, the Brillioun zone integration was 
performed using the Methfessel-Paxton scheme with a 
smearing of 0.2 eV and for the final structure, the energy 
was calculated using the linear tetrahedron method with 
Blochl corrections 

Our method of exploring the ternary phase space re- 
quires initial guesses of the crystal structure for a vari- 
ety of compositions. The results presented here take ad- 
vantage of the fact that a compound's structure is often 
reasonably close to that of a related ternary compound 
composed of different elements.'^ For example, the struc- 
ture of CeRhIn5 can be used as a guess for the structure 
of CeTMIus, where TM is any transition metal. Follow- 
ing this approach we assembled a set of candidate struc- 
tures from known compounds in the Ce-Ir-In systems 
and related binary and ternary systems. Specifically, we 
selected the known elemental, binary and ternary com- 
pounds of the Ce-Ir-In, Ce-Pt-In, Ce-Rh-In, and Ce-Pd- 
Al ternary and the Ru-In binary systems. We did not, 
however, include all known compounds of these systems, 
but rather selected those which experimental experience 
has shown most often form in Ce-Ir-In ternary phases 
and other chemically-similar systems. 

After determining the formation energy for all the can- 
didate structures, we determine the 3-dimensional convex 
hull, by plotting the composition versus formation energy 
with respect to the elemental phases. Those compounds 
that form the convex hull are energetically stable while 
the compounds above the hull have a positive formation 
energy against competing phases and are unstable and 
hence difficult or impossible to synthesize. From the hull, 
we can also identify competing phases. Calculating the 
formation energy of a compound with respect to its com- 
peting phases allows us to identify its degree of stability. 
Stable compounds will exhibit a negative formation en- 
ergy with respect to competing phases. 

The success of this project is due to the close communi- 
cation between experimental and theoretical efforts. The 
sharing of preliminary information between the groups 
about particular structures and entire portions of the 
phase diagram was instrumental in accelerating the dis- 
covery process. 
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FIG. 2: (color online) X-ray powder diffraction of Celr4ln. 
Tick marks indicate the expected peak positions based on the 
MgCu4Sn structure type. The difference between the mea- 
sured (red) and calculated (green) spectrum is shown above 
the measured spectrum in blue. The inset displays the crystal 
structure of Celr4ln (red spheres = Ce, blue spheres = Ir, and 
green spheres = In). 



C. Details of the Synthesis 

Compounds can be synthesized in a large variety of 
ways. For those compounds which were previously un- 
known to our knowledge (namely Celr4ln and Ce2lr2ln) 
we attempted to synthesize these crystals by arc melt- 
ing. The synthesis of Celr4ln was carried out via arc- 
melting the respective elements in a stoichiometric ratio 
of Ce:Ir:fn - 1:4:1. The resulting button was flipped over 
3 times and completely melted in order to ensure the re- 
action product was homogeneous. After the arc-melting, 
the pellet was weighed and was found to have less than 
a 5 mg loss. X-ray powder diffraction on the as cast pel- 
let was completed on an X-I Scintag diffractometer in 
the 20 range of 10 - 80°. The powder pattern was re- 
fined using the JADE software package. The attempted 
synthesis of Ce2lr2ln was carried out via arc-melting the 
respective elements in a stoichiometric ratio of Ce:Ir:In 
of 2:2:1. Characterization measurements of the magnetic 
susceptibility and heat capacity were performed in com- 
mercial systems from Quantum Design. 



structures, and the three elemental phases. The results 
are tabulated in Table [ij and the corresponding convex 
hull is presented in Fig. [T] Of the five known ternary 
compounds, three of them were computed to be stable. 
The two others, Ce2lrln8 and Celrln2, had only small 
positive formation energies of 0.7 and 40.2 meV with re- 
spect to competing phases, respectively. Of the eight 
unknown ternary structures, six have positive formation 
energies relative to competing phases, and thus we would 
expect them to be difficult to synthesize. However, two of 
these compounds have negative formation energies with 
respect to competing phases. Consequently we predict 
two new stable compounds: Ce2lr2ln and Celr4ln. Sub- 
sequent experiments described in Sec. |IIIB| successfully 
synthesized the Celr4ln. 

We characterize the electronic, magnetic and vibra- 
tional properties of the new Celr4ln compound. Phonon 
calculations for the cubic Celr4ln compound show an un- 
stable mode and predict a rotational distortion of the Ir 
tetrahedra. Relaxing the crystal structure along this dis- 
tortion mode rotates the tetrahedra by approximately 
11° and couples to a tetragonal distortion of the lattice 
parameters elongating the lattice parameter parallel to 
the axis of rotation. The space group symmetry is low- 
ered from cubic F43m (216) to tetragonal /4 (82). The 
lattice parameter change slightly from a = 7.566 A of the 
cubic structures to a = 5 = 7.508 A and c = 7.640 A. 
The symmetry breaking lowers the energy by about 300 
meV per unit cell. The X-ray diffraction pattern in Fig. [2] 
shows no evidence of the cubic cell undergoing any such 
distortion. This discrepancy between calculations and 
experiment may potentially be due to long-range effects 
not captured in the limited cell size of the calculations, 
R.g-i disorder among the Ce and Ir sublattices. 

Calculations of the electronic and magnetic structure 
of the distorted Celr4ln compound starting from vari- 
ous potential magnetically ordered states always result 
in a non-magnetic ground states. Each Ce atom con- 
tributes one /-electron to the bonding. The /-electrons is 
delocalized and no magnetic moments are present. This 
suggests that the magnetic features seen in the experi- 
ments described in Sec. [Ill B| are indeed a consequence of 
a small amount of impurity phases. 



B. Experiment 



III. RESULTS 
A. Theory 

For the Ce-Ir-In phase diagram we computed the 
formation energy of five known ternary structures 
(Ce5lr2ln4'i^, Celrln=|i^, CelrIrP, Ce2lrInJiS!^ and 
Celrlnji^, eight currently unknown ternary struc- 
tures (Ce2lr2ln, Celr4ln, Celr2ln4, Celr2ln2, Celrlns, 
Ce3lr4lni3, Celr2ln7, and Celr5ln2), nineteen binary 



Following the theoretical prediction of stable Ce2lr2ln 
and Celr4ln compounds, we attempted to synthesize 
these compounds by arc melting. The analysis of the 
X-ray powder diffraction data of Ce2lr2ln indicates mul- 
tiple competing phases which are not easily identifiable. 
Annealing of the Ce2lr2ln sample at 800° C for two weeks 
did not significantly alter the powder X-ray diffraction 
pattern. In contrast, the X-ray analysis of the arc melted 
Celr4ln sample shown in Fig. [2] confirms the formation 
of the predicted Celr4ln compound with lattice param- 
eters a = b = c = 7.579 A in good agrement with the 
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TABLE I: Calculated energies, competing phases, and space groups associated with Ce-lr-In phase diagram. The total energy 
of the compound, -Etotai, represents the calculated value per formula unit. The heat of formation of the compounds, A_Ef, 
relative to the elements and the formation energy relative to the competing phases, AiScomp, are given units of meV/atom. 
Compounds listed in bold are known to exist. 



No. 


Compound 
[eV/f.u.] 


-Etotai 

[meV/atom] 


[me V/ atom] 


^-£'comp 


Competing Phases 


Space Group 


1 


Ce5lr2ln4 


-64.081 


-588 


-65 


Cealus, Cesira, Cealrain 


Pbam (55) 


2 


Cealrain 


-35.881 


-749 


-42 


Ce5lr4, CeIra, Celrln 


PA/mhm (127) 


3 


Celrlns 


-30.347 


-395 


-35 


Cealrlug, CelrIn, Irlua 


P4:/mmm (123) 


4 


Celrin 


-19.415 


-689 


-32 


Cealrain, Celrlns 


P62m (189) 


5 


Celr4ln 


-46.651 


-455 


-8 


Cealrr, Irlna 


li (82) 


6 


Cealrlns 


-45.594 


-400 


1 


Celna, Celrlns 


P4/mmm (123) 


7 


Celrlna 


-21.987 


-520 


40 


Celrin, Celrlns 


Cmcm (63) 


8 


Celr2ln4 


-36.504 


-375 


44 


CeIra, Irlna, Celrlns 


Pmna (62) 


9 


Celralna 


-31.124 


-472 


56 


CeIra, Irlua, Celrlns 


P21/m (11) 


10 


Ce3lr4lni3 


-93.408 


-346 


70 


CeIra, Irlua, Celrlns 


Pm3m (223) 


11 


Celrlns 


-24.435 


-394 


89 


CelrIn, Celrlus 


PmSn (223) 


12 


Celrain? 


-43.520 


-197 


140 


Irlna, CeIra, Celrlns 


lA/mmm (139) 


13 


Celrslna 


-56.784 


-181 


189 


Cealrr, Ir, Irlna 


PA/mmm (123) 


14 


Celng 


-15.255 


-412 


-157 


In, Cealns 


Pm3m (221) 


15 


CeIra 


-25.994 


-782 


-149 


Ceslr4, CeIra 


FdSm (227) 


16 


Irlns 


-17.280 


-187 


-43 


In, Iralnr 


P42/mnm (134) 


17 


Ce5lr4 


-71.542 


-717 


-43 


Cesira, CeIra 


Pnma (62) 


18 


Cealns 


-33.643 


-382 


-30 


CeaIn, Celna 


Cmcm (63) 


19 


Cesirs 


-61.358 


-640 


-18 


Ceyira, Ce5lr4 


Pi/ncc (130) 


20 


Irlna 


-14.526 


-184 


-18 


Ir, Irlna 


Fddd (70) 


21 


Cesin 


-21.036 


-170 


-17 


Ce, Cealns 


PmSm (221) 


22 


Ceylra 


-73.297 


-520 


-7 


Ce, Cesira 


P6amc (186) 


23 


CeIra 


-34.913 


-602 


-5 


Cealrr, CeIra 


R3m (166) 


24 


Cealrr 


-78.684 


-535 





Ir, CeIra 


mm (166) 


25 


Iralnr 


-46.201 


-172 


13 


Irlua, Irlna 


ImSm (229) 


26 


\-jt:a J.1 


—98 959 


-399 


34 




Pnmn (f\9') 


27 


Ceain 


-14.929 


-168 


49 


CeaIn, Cealns 


P63/mmc (194) 


28 


Ceirs 


-52.303 


-347 


54 


CeaIrT, Ir 


PQ/mmm (191) 


29 


Celna 


-12.050 


-334 


58 


Cealns, Celna 


Imma (74) 


30 


Irin 


-11.532 


-58 


80 


Ir, Irlna 


C2/m (12) 


31 


Irain 


-28.895 


59 


128 


Ir, Irlna 


Pm3m (221) 


32 


Iralna 


-31.479 


42 


152 


Ir, Irlna 


P31c (163) 




Ce 


-5.933 








FmSm (225) 




Ir 


-8.857 








Fm3m (225) 




In 


-2.558 








M/mmm (139) 



calculated lattice parameters. 

The crystal structure of Celr4ln is presented in the in- 
set of Fig. [2| The Ce sublatt ice exhibits an fee structure 
composed of corner sharing tetrahedra, which satisfies 
the condition for geometric frustration of the magnetic 
moments on the Ce atoms. Consequently, we measured 
the magnetic susceptibility and heat capacity to deter- 
mine the basic properties and whether there is evidence 
for magnetic frustration in Celr4ln. 

The magnetic susceptibility is shown in Fig. [3j An in- 
dication for ferromagnetic ordering at low temperature 
is indicated by the rapid rise in the susceptibility. The 
inset showing M vs. H further supports ferromagnetic 
ordering at low temperature albeit with a small ordered 
moment. The inverse susceptibility shows a large region 
of curvature up to room temperature. Similar temper- 
ature dependence seen in SmPt4ln was interpreted in 



terms of a Van Vleck contribution plus crystal electric 
field (CEF) effects.'^Sl However, the CEF ground state for 
Ce in Celr4ln is most likely magnetic. Alternatively, the 
large region of curvature could indicate a relatively large 
energy scale for the onset of exchange interactions, for 
which the ordering between moments is suppressed to 
much lower temperatures. 

The heat capacity data shown in the inset of Fig. [3] 
shows a clear albeit broad magnetic ordering transition 
at 2.5 K consistent with the susceptibility data and resis- 
tivity measurements (not shown) and a smaller anomaly 
at 6 K. The field dependence of the anomaly at 2.5 K 
(not shown) and the shape of the transition in suscepti- 
bility again suggests that the magnetic ordering is ferro- 
magnetic. The extrapolated zero temperature Sommer- 
feld coefficient in the ordered state is 55 mJ/(mol-K^). 
From the specific heat data, we calculate the entropy 
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FIG. 3: (color online) Magnetic susceptibility of Celr4ln ver- 
sus temperature. Upper left inset displays the specific heat of 
Celr4ln versus temperature and the lower right inset shows 
the magnetization as a function of magnetic field at 2 K. 



S ^ J C/TdT. However, only about 5% of i? In 2 in en- 
tropy is recovered at the transition temperature, much 
smaller than expected for a system of local magnetic mo- 
ments at the ferromagnetic ordering temperature. 

The small amount of entropy and small magnetic mo- 
ment inferred from the susceptibility measurements sug- 
gests that these features are a consequence of a small 
amount of impurity phases. We note that there are small 
unaccounted peaks in the X-ray pattern shown in Fig. 2. 
While we believe this is the most likely explanation a 
few alternative scenarios still exist. First, the reduced 
entropy may be a consequence of Ce moments having 
largely itinerant character {i.e. in a Ce^+ configuration) 
as has been suggested to occur in CeNi4ln.^ Alterna- 
tively, a large majority of the magnetic entropy may 
be quenched at much higher temperatures as a result 
of disorder-, RKKY-, and/or geometric-induced mag- 
netic frustration effects as has been su gge sted to occur in 
the i?Cu4ln {R = Gd-Er) compounds. Finally, Kondo 
screening is another alternative to quenching the mag- 
netic entropy at higher temperatures without long range 
order. 

Of all the CcTM4ln compounds, CePt4ln is most sim- 
ilar in structure and chemistry to Celr4ln. CePt4ln is 
a local moment system displaying heavy fermion behav- 
io] |23|24 |^ with a Kondo temperature of ~10 K and short 
range ordering occurring at 250 mpp^l which would favor 
some combination of the latter two scenarios for explain- 
ing the reduced entropy in Celr4ln. If Ir hybridizes more 
strongly than Ft, this would favor the scenario where 
the magnetic anomalies are a result of impurity phases, 
and Celr4ln is approaching mixed valent behavior akin to 
CeNi4ln. Consequently, it is currently inconclusive as to 
whether geometric frustration is playing a significant role 
in Celr4ln. More measurements are needed to determine 
the origin of the physical properties in Celr4ln. 



IV. DISCUSSION 



In retrospect, it may not be surprising to anticipate 
the discovery of Celr4ln when the isostructural analog 
CePt4ln is known. However, all of the unknown Ce-Ir-In 
compounds computed in Table [l] contain closely related 
known isostructural analogs. Furthermore, it is remark- 
able that these relatively simple DFT calculations which 
neglect zero point motion, finite temperature, and elec- 
tronic correlation effects, accurately predict the existence 
and lack thereof, for 10 out of the 13 ternary compounds 
and 15 of the 19 binary compounds investigated here 
in the Ce-Ir-In ternary phase diagram. In addition, for 
Ce2lrln8 and Celrln2 which are predicted to be unstable 
but are known, and Ce2lr2ln which is predicted to be 
stable but has yet to be found, the formation energies 
relative to competing phases are sufficiently small that 
we might expect that some of the neglected effects men- 
tioned above might change the conclusion with regards 
to compound formation. We must then also acknowl- 
edge that given the small positive formation energy with 
respect to competing phases for Celr2ln4 and Celr2ln2, 
which is less than 56 meV/atom, it would not be surpris- 
ing if either or both of these compounds are eventually 
synthesized. 

In addition to the possible shortcomings of the DFT 
calculations mentioned earlier, we also recognize that the 
phase space of potential compounds which we are ex- 
ploring is remarkably small. Not only have we selected 
a relatively small number of compositions, but we have 
also assumed a particular structure type for each of these 
compositions. In this respect, the addition of simulated 
annealing, genetic algorithms and other approache^^H to 
this approach may aide in discovering additional struc- 
tures not considered based on analogy. 

The fact that arc melting Ce2lr2ln and annealing did 
not produce the anticipated structure does not neces- 
sarily mean that the compound is unstable. However, 
it shows that Ce2lr2ln can not be a congruently melt- 
ing compound and that the phase is either not stable at 
800° C or that a direct transformation from the compet- 
ing phases may exhibit slow kinetics. Alternative synthe- 
sis routes should be considered to stabilize the compound. 
These could be as simple as using a different temperature 
annealing profile or using another synthesis route such as 
fiux growth.^ Thus, we anticipate that the effort in try- 
ing to synthesize a theoretically predicted compound by 
this method will be commensurate with the interest that 
the expected properties of that compound would gener- 
ate. 

Despite these potential complications, the fact that 
our combined computational and experimental approach 
demonstrates an accuracy of better than 80% in identi- 
fying which compounds exist or not, proves that this can 
be a useful method for computationally identifying new 
stable structures. 
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V. SUMMARY 

We present a simple method combining density- 
functional calculations with arc-melting synthesis to effi- 
ciently investigate the structural stability of hypothetical 
compounds with the aim to improve the search for new 
compounds with desired structures and properties. As a 
proof of principle we applied our approach to the Ce-Ir- 
In phase diagram and predicted two previously unknown 
compounds, Ce2lr2ln and Celr4ln, to be stable. Sub- 
sequent attempts to synthesize these compounds led to 
the discovery of Celr4ln. We expect that this approach 
will lead to the discovery of many more interesting com- 
pounds with specifically desired properties. 
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